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1. Introduction 



For three decades our microscopic theory-QCD, has been unable to describe soft (long distance) interac- 
tions. It happens, partly, because of the embryonic stage of our understanding of the most challenging 
problem in QCD: confinement of quarks and gluons. We believe, however, that the main features of high 
energy scattering at long distances can be described in the framework of QCD, and such a description is 
based on the phenomenon of saturation of the parton densities in QCD |3|. Indeed, at short distances 
the parton densities increase with the growth of energy due to gluon emission, which is described by the 
QCD linear evolution equations ||. Such an emission is proportional to the density of emitters (partons). 
The process of annihilation is suppressed since it's strength is proportional to the parton density squared, 
reflecting the fact that two partons have to meet each other in the small volume of interaction. However, 
at high energy the density of partons becomes so large that the annihilation term tends to be of the same 
order as the emission term. It leads to a slow down in the increase of the parton density, and the system 
of partons lives in the dynamical equilibrium with the critical value of density (it is saturated). We have 
several theoretical approaches for a QCD description of the parton system in the saturation domain (see 
Ref. [Q, ^ 1^, ^, 1^, 1^, ^, |ll|, but the equations that have been proposed turn out to be so complicated, 
that our knowledge about this domain mostly stems from numerical solutions of these equations p 



The alternative approach that has been developed during the past decade is to build models [0, |15|] that 
incorporate the main qualitative properties of these solutions, such as the geometrical scaling behavior of 
the scattering amplitude ||l9| , 20, ^] and the existence of the new dimensional scale (the so called saturation 



scale) ^, ^ , which increases in the region of low x (high energy) . Such models successfully described all 
the data on deep inelastic scattering, including the behavior of inclusive and diffractive cross sections, at 



a low value of photon virtuality |14]. 



The main goal of this paper is to describe the processes that occur at long distances (soft processes) , 
where we cannot apply perturbative QCD, using the saturation model. In this model all the needed 
phenomenological parameters, in particular the energy behavior and the value of the saturation scale, were 
fitted using the DIS data. The main ingredient of the high energy phenomenology during the past four 
decades, namely the soft Pomeron, is not introduced in this model. In simple words, we replace the soft 
Pomeron by the scattering amplitude that describes the behavior of the saturated parton system. The idea, 
that the matching between soft and hard interaction occurs in the saturation region, and can be described 
by the amplitude that originates from the high density QCD approach, is not new. This idea has been in 



the air for a number of years and have been advocated in several lectures on the subject [22|. The first 
attempt of a practical application of this idea was done in Ref. |2^ (see also Ref. [^]). The result of 
this paper was encouraging since the estimates gave a reasonably good agreement with the experimental 
data on the total cross section of the pion, kaon, and proton interaction with the proton target. In this 
paper we continue to develop these ideas and we give the comparison with the experiment data for the full 
set of the experimental data of soft interactions, that include the energy behavior of the total, elastic and 
diffraction cross sections as well as a t- behavior of the elastic cross sections, and the mass behavior for the 
cross sections of diffraction production. 
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2. The main idea 



As has been mentioned, our main idea is to describe the so called 'soft' interaction in the high parton 
density QCD (hdQCD). In hdQCD we are deahng with the dense system of partons which cannot be 
treated in the perturbative QCD approach. Being non-perturbative in its nature, hdQCD leads to typical 
distances r = 1/Qg(x) <^ R ^ 1/Aqcd where R is the size of a hadron. The saturation momentum Qsix) 
is a new dimensionful scale which increases with the growth of the energy ( or at X ^ 0)|l|, I, |]. Therefore, 
the QCD coupling as {Qs{x)) is small and, because of this, we are able to develop a theoretical approach 
for such a dense system having only a limited input from the unknown confinement region. To illustrate the 
limitation of our approach, it is enough to look at the expression for the saturation scale in hdQCD [||, |^ 

Q'iiY) = Q^^(Y,)exp(as^^iY-Yo) - , ln(y/yo: 

(2.1) 

In Eq. ( ^ ) the energy dependence of QsiY = ln(l/x)) is predicted from hdQCD, (x(7) is the Mellin 
transform of the BFKL kernel) but the scale , Q'^{Yq) stems from the confinement region, demonstrating a 
need for some input from this theoretically unknown region. Our hope that such input will be a limited 
number of constants. 

Let us consider the deep inelastic electron-proton scattering (DIS), to illustrate our point. If the 
virtuality of the photon is very large, ^ Qg(x), the typical distances r ~ l/Q, and we can safely 
apply the perturbative QCD approach based on the operator product expansion, and the DGLAP evolution 
equation. However, if Ajg^jjj <^ <^ Q'^{x), the situation changes crucially. All terms in the operator 
product expansion become of the same order, and we have to use the hdQCD approach. In this approach, 
the cross section a^*P shows a geometrical scaling behaviour 1 27], namely, a'^*^ = (^1/Ql{x)^ F ((5^/Qs(x)) 




which has been confirmed experimentally. This means that the typical distances r ~ 1/Qs{x) <C 1/Aqcd, 
which are short. 

Our idea is that we can describe « Ag^^,, in the framework of the same hdQCD approach. 
Having the geometrical scaling behavior in mind, such an idea does not look crazy. However, in DIS we 
are dealing with the total cross sections which are related to the amplitude integrated over the impact 
parameters (6). For a treatment of the soft interaction observables, we need to know the b dependence. 
As was shown in Ref. [^8[, hdQCD predicts the power-like decrease at large value of 6, which contradicts 
both the theoretical estimates and experimental observations. Therefore, we need an input from the 
confinement region to specify the b dependence of the scattering amplitude. Let us first discuss large 
b > 1/Aqcd 3> l/Qs{x). The scattering amplitude for DIS A (^Q'^,W;t) depends on the photon virtuality 
Q^, energy W = y/s = Q/ y/x and on the momentum transferred t = —q^ < 0. Using the dispersion relation 
in the t-channel we have 

.2 w..^ _ 1 r dt' ImtA{Q^,W;t') 



A{Q^W;t) = - / Vl"; ' (2-2) 
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Calculating the amplitude in the impact parameter representation we can reduce Eq. ( |2.2| ) to the form 
A{Q\W-h) = I d'qA{Q\W;t)e^'^-^' = ^ j^^ dt' Im,A {Q^W;t') J ^^^^^ (2.3) 

= ^ dt'Ko(^Vt^bj ImtA{Q^,W;t') ^^^Zi^T^) C ImtA {Q\W;t' = Ami) Ko{2m^b) 

where in the constant C, all the numerical factors have been absorbed. 



In Eq. (|2.3| ), the 5-dependence is determined by the mass of the lightest hadron (pion), and cannot 
not be reproduced in perturbative or/and high density QCD. However, one can see from Eq. ( |2.3D that we 
have a factorization of the b dependence, and the dependence on the photon virtuality and energy. 

As it is well known (see Ref. fl^ and references therein), in the kinematic region as < as ln(l/x) < 
1/as in the saturation region the hdQCD approach reduces to the BFKL Pomeron calculus. We generalize 
Eq. ( |2.3D for the BFKL Pomeron in the following way 

short distances long distances 

Abfkl {Q^,oc;b) = A BFKL {Q^ ,x;t = (S) X S{b) (2.4) 



One can see that Eq. (|2.4| ) claims a kind of factorization between short and long distances. The short 
distance part we can describe in hdQCD, while the long distance part has to be taken from the confinement 
domain. In a more general way, we can re-write Eq. (^^) in the form 

short distances long distances 
ABFKL{Q\x-b) = I d^b' ABFKL{Q^x;b') X S(6-6') (2.5) 

where b' in the short distance part is of the order of b' ~ 1/Qs{x). 



The approach of Eq. ( |2.4| ) and Eq. (|2.5D , at first sight contradicts the high energy phenomenology based 
on the soft Pomeron and Reggeons, as well as well as the lattice calculations [^]. In both approaches, 
the Pomeron trajectory has a'jp > 0, and at positive t the glueballs lie on this trajectory. In Eq. ([2.4]) 
and Eq. ( |2.5| ), the Pomeron slope a'p = 0. We will show below that we will be able reproduce the 
experimental data on energy dependence of the elastic slope which used to consider as the argument for 
a'jp = 0.25 GeV~^. On the theoretical side, the only theory which can treat at the moment on the same 
footing both short and long distances: N=4 SYM with AdS/CFT correspondence, leads to a picture which 



is in striking agreement with our approach [ 29 1. Namely, in this approach at t > ( resonance region) we 



have a normal soft Pomeron with a'jp > while at t < (scattering region) a'jp = 0. 

3. The model 
3.1 Motivation 

Our building of the model is based on the dipole approach to high energy QCD [^]. In this approach, 
the evolution of the parton (colorless dipole) system can be written in the form of the equation for the 



- 4 - 



generating functional, with a transparent probabilistic interpretation for them. The generating functional 
is defined as [31| 

/n 
Pn{Y - Yo; xi,yi;...;xi,yi;...;xn,yn) y^) t/i (3.1) 



n=l i=l 



where u{xi,yi) = Ui is an arbitrary function of Xj and y^. The coordinates {xi,yi) describe the colorless 
pair of gluons or a dipole. Pn is a probability density to find n dipoles with the size Xi — yi, and with 
impact parameter (xj + yi)/2. Assuming that we have only a decay of one dipole to two dipoles, directly 



from the physical meaning of P„ it follows |31, 32 



d PniY] . . . ] Xi, yi] . . . 'i Xji, y-n) (ri n\ 

dY " ^"^-^^ 

= ^ Vi^2<^{Pn-i{Y;...]Xi,yi;...;Xn,yn) - Pn{Y;...;Xi,yi;...;Xn,yn)) 

i 

where denotes all necessary integrations and Vi_>2 is the vertex for the decay of one dipole to two 
dipoles. This vertex is equal to 

yi_2((x,y) ^ (x,z) + (z,y)) = J^)2(^^_^)2 (3-3) 

Eq. ( |3.2| ) is a typical Markov's chain which takes into account the s - channel unitarity on each step of 
evolution since it has two terms: the birth of the new dipole due to the decay of the parent dipole (positive 
term in Eq. ( p. 2D ) and the death of the dipole due to the same process. 



Eq. ( p.2| ) can be rewritten as the equation for the generating functional |32, namely, 
dZ {Y-Yo;[u]) _ /'^2^^2^^2^^^^^((^^^)^(^^^)^(^^^)) d ziY-Yo;[u]) (3.4) 



dY J du{x,y) 

Z (Y — Yq;[u]) satisfies the initial and boundary conditions: 

Initial conditions: Z {Y — Yq = 0; [u]) = u{x,y) ; (3-5) 

Boundary conditions: Z {Y — Yq; [u = 1]) = 1 ; (3-6) 

The advantage of the generating functional is that we can write it in terms of this simple functional 
formulae for the scattering amplitude, which has an obvious partonic interpretation. Indeed, the scattering 
amplitude in the lab. frame can be written in the form ||6|, |3^, 33 1 

N{Y-Yo;r,h) = (3.7) 
^ ^, / (f yi Pn (r, h; {x^, yj, Y - Yq) 7„ (r, b; {xj, yj, Yq) 

n=l i=l 
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where Y — Yq » 1 but Yq ^ 1 and 

n ^ 

Pn{Y -Yo,{xi,yi}) = YIyI (^-^o;{w(a;i,yi)})U(x„j/0=i (3-^) 

ou[Xi,yi) 

and 7„ is the amphtude of interaction of n dipoles at low energy ( small values of rapidity Yq) with the 
coordinates Xi,yi with the target with size r and impact parameter b. In the parton language the generating 
functional determines the parton wave function for which we have the evolution equation Eq. ( |3.4| ) in QCD. 
7„ can have a non-perturbative origin. We assume that 7„ (r, 6; {xj, yj}, Yq) = IliLi -^o)- This 

assumption means that the dipoles inside the target have no correlation, which is correct, as far as we 
know, only for a nucleus target. 

In spite of the transparent physics behind our equations, the equation for the generating functional 
is difficult to solve analytically, especially if we generalize Eq. ( |3.4D to the case of the so called Pomeron 
loops ^ |3^. As has been mentioned, we wish to find a model solution to Eq. (|3.4|) . We will do 



this using two key observations of Ref.[12|. The first one is the fact that in the huge kinematic region of 
asY < 1/as, the system of partons that we are dealing with turns out to be a system of non- interactive 
partons. In other words, it means that the generating functional of Eq. ( |3.1| ) can be rewritten in a simpler 
form, namely. 

Z {Y -Yo;{u{xi,yi)}) = ( / d^x' (fy' P{Y -Yo;x,y;x' ,y') {u{x' ,y') - l) ] (3.9) 



^ ^ / d'x'd'y'P{Y-Yo;x,y;x',y') {u{x',y') - l) (3.' 

n=l ^ 



where P{Y — Yq; x, y; x' , y') is the amplitude for one BFKL Pomeron exchange, between dipoles {x, y) and 
{x',y'). 

In Ref.[12| it was shown , using the analytical solution of Ref. [34| , that actually C„ = 1 for the 
dipoles with the size r ~ l/Q^ where Qs is the saturation momentum. 

Therefore, Eq. ( |3.9| ) can be simplified and rewritten in the form 

Z iY-Yo;{u{xi,yi)}) = exj> (^J d^x' d^y' PiY - Yo; x,y; x' ,y') {u{x' ,y') - l)^ (3.10) 

and the amplitude has the form 

N{Y-Yo;x,y) = 1 - exp j d^x' d^y' P {Y - Yq; x,y; x' ,y') ^ {Yq; x' ,y')^ (3.11) 
where 7 {Yq; x', y') is the scattering amplitude of the dipole {x\ y'), with the target at low energy. 
3.2 Main formulae and assumptions 



Eq. ( 3.11 ) is the main formula for constructing our model. The first ingredient of the model is the choice 



of the amplitude 7 and the value of Yq. We choose this amplitude in the form [^] 

^{YQ;x',y') = ^mo4Yo;x',y') = ^^^^^ xqG^^^^^ {xq, p') S{b) (3.12) 
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where r = x — y and h = ^{x + y) with xq = 10 ^. xG in Eq. (3.12) is the solution of the DGLAP 



evolution equation which describes the DIS experimental data, /i is equal to 

C 



fi' = 4 + ^ (3.13) 



r 



where /xq cind C are phenomenological parameters which has to be found fitting the experimental data 
on DIS (see Ref.[15| for details). Eq. ( p.l3| ) means that "ymod {Yq] x' ,y') at large dipole size rfiQ 3> 1 is 



determined by the following expression 



lmo4Yo;x',y') = ^^l^^r'xoG^^^^^ixo,^l)Sib) (3.14) 



Of course Eq. ( 3.14| ) has no theoretical basis, and our hope is that the amplitude at low x will be not very 



sensitive to this kinematic region, in the initial condition. 

S{b) is the impact parameter profile function, which is chosen as the Fourier image of the proton form 
factor, namely, 

where R is the proton radius {R = 0.89 /m). 



In Eq. (3.11) P {Y — Yo;x,y;x' ,y') stands for the BFKL Pomeron Green's function. This Green's 
function describes the energy (x) evolution of the gluon system in the region of high energy (low x). This 
evolution allows us to find the dipole scattering amplitude or parton density at low x from the initial 
condition: the amplitude at lower x = xq, but for any value of the dipole size (r). 

However, we know that the BFKL Pomeron alone cannot describe the experimental data, since it 
determines the anomalous dimension only at low x. The experimental data shows a good agreement with 
the DGLAP evolution equation, which is written as the evolution in ln(//^) (see Eq. ( 3. 12] )), which gives 



us the parton density from the boundary condition: the parton density at = fi^ but at any value of x 
including the region of low x. Strictly speaking, it cannot be used in Eq. ( 3.11| ), but assuming Eq. ( p. 131 ) 



for the scale /U^, we impose the condition that at long distances, the typical momentum scale in the linear 
evolution equation is equal to fi^. In the framework of this assumption, it looks reasonable to assume that 



n{Y,r,b) = 2 j (fx'd^y' P{Y -Yo;x,y;x',y')j{Yo;x',y') 

= ^llMlll xG^«^-^^(x, /x) S{b) (3.16) 
3 

Finally, we can write the following model expression for the scattering amplitude; 

N{x;r,b) = 1 - exp f -^17 (x, r, 6) j (3.17) 



where (x, r, b) is given by Eq. ( 3.16| ). 
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3.3 Fixing the values of the phenomenological parameters 



In Eq. ( 3.13| ) we have introduced two phenomenological parameters: fiQ and C. The first determines the 
value of the virtuality of the probe in DIS processes from which we can start using the perturbative QCD 
evolution equations. The second (C) gives the relation between scales in the momentum representation 
and in the coordinate representation. In the case of the DGLAP evolution, C = 4 ( see Ref. |38]), but 
since the amplitude has a more complicated form than Eq. ( |3.17 ), the value of C can be different. 

Two more parameters stem from the initial condition for the gluon density in Eq. ( 3.161 ), namely, 



xG 



DGLAP 



") Qo) 



A 



{l-xf 



(3.18) 



which corresponds to the exchange of the soft Pomeron with the intercept uiq at the initial hard scale 
jx^ = Qo traditional high energy phenomenology, and the factor (1 — x)^ reflects the valence quark 
dependance. In our approach we do not assume that the soft Pomeron exists, and we view Eq. ( 3.18| ) just 
as the simplest function that reflects the behavior of the experimental data in DIS. 



We use the observation of Ref. |39], that the anomalous dimension in the leading order can be written 
in a very simple form, namely, 

7(c^) = asi- - l] (3.19) 

and with the initial 



(3.20) 



The explicit solution to the DGLAP evolution equation with 7(w) of Eq. ( p.l 
condition of Eq. ( 3.1^ ), has been found in Ref. [15] and looks as follows 



xG 



DGLAP, 



fc=0 ^ ^ 



where 



!;G(^')(y = ln(l/x),t,cufc) = Ae'^ + '^^y 



dy' e 



-t^k y 



/i 20?7 + 1 



with 



Uk = Wo 



t 



4jVc In (/^VA^) 
6o ^n(/.g/A2) 



2nf 

11 and y = ln(l/a;) 



(3.21) 



(3.22) 



In Eq. ( 3.16| ) the typical scale of hardness in the gluon structure function is determined by the process 
of gluon emission while the factor as r'^ takes into account the integration over the wave function of the 
dipole with the size r. Having this in mind, we introduce a different scale for as in Eq. (3.16). Finally, we 
use 

7r2 as{jl^) 



r2(a;, r, h) 



■r^xG^^^^^{x,ij)S{b) with fl 



,~,2 



~2 C 
^0 + ^2 



(3.23) 



where xG^'-^^^^{x, ^) is given by Eq. (|3.21|) and the initial hard scale /Xg is different from that of fi^. We 
perform a fit of the saturation model, to the DIS experimental data and fixing the parameters, we wish to 
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Model 


4 


4 


Ql 


C 


Wo 


A 


xVd.o.f. 


A 


0.23 


1.23 


1.0 


1.0 


0.028 


1.81 


1.04 


B 


0.069 


1.23 


0.58 


1.0 


0.087 


0.6 


3.08 



Table 1: Resulting parameters from the fit to DIS data Model A gives a very good description of the DIS data, while 
model B reproduces quite well the data on soft processes. 



describe the soft data without using the concept of the "soft" pomeron. The fit procedure was performed 
by using the minuit routine, the statistical and systematical errors have been added in quadrature. We 
have used all of the recent data on deep inelastic scattering processes, from different collaborations. The 
most suitable were HI [16| and ZEUS |17, It was observed, that HI data for medium and large values 
of coincides with that from the ZEUS col., but up to a scaling factor of 1.05. Hence, for the fitting 
procedure, we took only the ZEUS data for vary small, medium and large values of Q^. Since we are 
interested in the high energy description, we took data below x < 0.01 and 0.045 < < l50GeV^. The 
upper limit on virtuality is originated from the upper limit on the x variable. Fitting with HI scaled 
data, yields almost the same result. The total number of experimental points was 170. During the fitting 
procedure, we observed a strong correlation between different parameters, so we decided to fix a number 
of them. We choose to fix a parameter C and the initial hard scale Qq. The masses of the quarks were 
taken as follows: mught = 0.25 GeV and rricharm = 1-3 GeV. this parametrization is ~ 1.05. The 

summary of the fitting procedure is given in the table H 

The bolded font corresponds to the fixed parameters C and Qq, which were chooses to be equal to 1.0. 
The value of the hard scale ^Ug, ptg. is in the units of energy GeV. 

As one can see, x^/d.o.f. is close to 1 and fit gives very good description of all experimental data on 
DIS (see Fig. |l]). However, we present in Table 1 a second fit which leads to worse /d.o.f. but reproduces 
the data on soft interaction quite well as we will see below. It should be stressed that both fits describe 
the DIS data at large values of photon virtualities with small x^/d.o.f. . Therefore, we can consider model 
the A as an attempt to describe the experimental data on DIS and soft processes, assuming that Eq. ( ^l7|) 
is correct in the saturation region. In model B we explore a different idea: Eq. ( |3.17 ) is only approximate 
formula that incorporates the main qualitative properties of scattering amplitude for both DIS and soft 
scattering in the saturation region. The question which we try to answer using model B is the following: 
is it possible to give the unique description of long distance physics both in DIS and soft interaction based 
on the QCD amplitudes at short distances. 



4. Cross sections of hadron-hadron interaction at high energy 
4.1 General approach 



Eq. ( 3.17 ) gives a smooth continuation to the long distance physics describing the DIS data at very low 
values of the photon virtualities. We wish to extend this description to the hadron-hadron interaction 
without assuming something in addition, for example, the existence of the soft Pomeron. The simplest 
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Figure 1: Description of the DIS experimental data in our models. Dotted line describes the model B while solid 
line corresponds to Model A. 



formula we can write for the total cross section, is a straightforward generalization of the formula for the 
DIS cross section with the replacement of "^-y*{r) — > hadroniiri}) , namely 

/„ n n 
dFh / J] d\i \'^hadron{{ri})? N{x,rf,b) (4.1) 

i=l i=l 

where n is the number of dipoles that we need to introduce to describe a hadron. For example, for a meson 
we need only one dipole, while for the proton we have to introduce at least two colorless dipoles. 

The total elastic cross section can be written in the form: 

(TL Tl \ ^ 

/ n ^'^i E \'^hadron{{ri})\^ N{x,ri;b)] (4.2) 
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For the differential elastic cross section we have the following expression 



dad 
dt 



1 



(2vr)^ 



„ n n N 

Jo{qh)dH / n d\i 1^ 

hadron 



(4.3) 



i=l i=l 

where Jo is the Bessel function. It is even easier to calculate the slope in t at i = 0. It is equal to 



el 



d In da el / dt , 
df 



\t=o 



1 J b^rf^b /nr=i d^r, Zti \'^hadroni{r^})\' N{x,r,;b) 

2 I d^b fU-=id^nY:^=, l^hadroniinW N{x,n;b) 



(4.4) 



The process of diffractive dissociation is a more complicated phenomenon. Indeed, the eikonal type formula 
of Eq. ( p. 17 ) leads to diffractive dissociation in the state of n free dipoles which, being a system with a 
limited number of dipoles, cannot create a system of hadrons with large mass. In other words, this 
diffractive production falls down as a function of produced mass. Therefore, we can calculate diffraction 



in the region of the small mass using the Good- Walker formula [ 41 1 , namely 



low M 



(4.5) 



hadroniin])?' N{x,ri]b) 



/.. n n ^ I „ n n 

d% / n d^n \'^hadron{{ri})? iV'(x,r,;6) - / d'M / n ^'^^ E 
i=\ i=l i=l i=l 

the first term is the cross section for the production of the system of n free dipoles, while the second is the 
elastic cross section. We subtracted this term to find the cross section of the hadron state, which is different 
from the initial one. . For large mass diffraction, we have to develop a new approach. The simple form of 




Y 



Y 
Y{ 

y; 




Figure 2: Simple Pomcron diagrams for the total cross section (Fig. ^-a) and for the diffractive production in the 
region of large mass (Y ~Yi ^ h\{M'^ / sq) , Fig. |^b). The dashed line shows the cut Pomeron. 



the generating functional of Eq. ( 3.10 ) stems from the fact that the system of interacting Pomerons can 
be reduced to the exchange of non-interacting Pomerons after integration over rapidity Yi, as it is shown 
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in Fig. ||-a. In the case of difFractive production we can also calculate the Pomeron diagrams with one cut 
Pomeron (shown by dashed line in Fig. ||-b) but the value of y — Yi is fixed, namely, Y — Yi = ln(M^/so) 
where M is the mass of the diffractively produced system and sq is the energy scale (sq ~ IGeV'^). In more 
complicated diagrams, we have to integrate over rapidities (for example, over Y( and Y^ in Fig. |-b). For 
the Pomeron with the intercept larger than 1, as in our case, these integrations result in three diagrams in 
Fig. |2|-b. It should be stressed that the integral over the produced mass gives the dominant contribution 
at y — Yi = ln(M^/so) ~ ^/ots- However, it has not been taken into account in our eikonal - type model, 
which only describes the G-W contribution to the region of small mass. 

Therefore, for da/dM"^ at large ^ sq, for one dipole {x,y) we can write a generic formula (see 
Refs. IS 



M 



da 



high M 
diff 



dM^ 

l-exp(-n{M^;x,y;x',y') 



V,^2{{x',y')^{x',z) + {z,y')) 

1 



exp 



exp(-0(s;x,y)) 
{n{s/M^;x',z) +n{s/M^;z,y')} 



(4.6) 



exp ( --n{s;x',y') 



where denotes all needed integrations, and is a new opacity for the interaction of two dipoles: (x, y) 
and (x',y'), which we can build using the same approach as in Eq. (13.161) . We start to analyze Eq. (gj) 
considering the emission of an extra gluon (see Fig. ^). 




Figure 3: The process of diffractive production in 
the region of large mass in perturbative QCD. 




Figure 4: Feedback to the process of elastic scat- 
tering of dipole(x, y) due to the emission of an extra 
gluon. 



One can see that our process has three well separated stages. The first one is a penetration of the dipole 
(x, y) through the target without inelastic interaction. We introduce the factor "exp (—5^)" to describe this 
stage. This factor sums all Pomeron diagrams with Pomerons that carry the total energy of the process( 
rapidity Y). See the second diagram in Fig. |2|-b. The second stage is the emission of one extra gluon. The 



dipole decay is responsible for this stage with the probability given by Eq. (3^). The last, third stage is 
the interaction of two produced dipoles {x',z) and z,y') with the target. This stage is taken into account 
in Eq. ( [4.6| ) by the factor in the brackets. This factor has been discussed in Ref. |43|. 



In the first diagrams (see Fig. |^ the second stage is very simple: it is just the perturbative emission 
of one extra gluon. The simplified formula for this diagram has the same structure as Eq. ( [4. 61 ), but with 
a simple expression for exp (— r2(s; x, y)) i 1 — exp ( — fi(M^; x, y; x', y'))), namely 
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^ ^ 

2 



exp ( - <{ 17(s/M2;r;6) - ^^^(s/M^ n; 6) 



ro (1 - exp(-f](s/so;ri,6)))^[> (4.7) 



27r (r*! — r)^ 

where n is the number of dipoles in a proton and r = x — z . In Eq. ( [4.7D we consider the region of 
integration where ri <^ r. The second term in curly brackets takes into account a change for elastic 
scattering of the dipole ri due to the emission of one extra gluon (see Fig. ^.) Eq. (|4.7| ) is written in the 
leading log approximation of perturbative QCD in which we consider Q;5ln(M^/so) ~ 1 while as <C 1. 
Factor (l — exp (— r2(s/M^; r) + ^Q,{s / M"^ , ri; b))) is the amplitude of gluon-dipole scattering that has been 



discussed in Ref. The formulae for diffraction production is well known (see Refs. [44| for details). In 
the same approximation Eq. ( [4.6[ ) has the form 

^^high M 



= ^ d% n{d\, \^hadron{{ri}t {e^Y>{-^{s/so;n,b))) (4., 

•' i=l 

I d\r^ jd%' (l-exp(-0™(MV.o;r„r))) j^^ (l - exp (-1](s/M2; /; 6)))', 



where Q^^^^ (M^/so;rj,r) is the solution of the BFKL equation with the initial condition: 



This O is equal to 



dn-^^^{Myso = l;r,,r) _ a. . 



Q^P^HM^/s^ = l-r,,r) = / ^2_e-WWMV.o)+7in{r?/r^) (4.10) 
^ ^ J 27riuj{j) 

r2 a;(0) Y a;"(0)ln(M2/so) 

where the eigenvalues of the BFKL equation uj{iy) can be found in Ref. 

The factor 1 — exp (—0"^^^^) describes the inelastic cross section which in terms of the Pomeron 
diagrams of Fig. ^, reflects the possibility that in the fourth of Fig. ^b, we can have two cut upper 
Pomerons. 

We would like to stress again that Eq. ( 4.10| ) describes the dependence of the diffractive production 
cross section, in the region of large mass, while Eq. ( [4.5| ) is written for the low mass diffractive cross section. 
Therefore, to calculate the cross section for diffractive production, we need to calculate 

^d^ff = cj'rff'HEq. (P)) + a'^ff'HEq. (H)) (4.11) 
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4.2 Hadronic wave functions 



As seen from the formulae in tiie previous subsection, we need to know the hadronic wave functions. We 
have made an assumption in writing these formulae, that the correct degrees of freedom at long distances 



are the colorless dipoles at least at high energy[^, 24|. This is one of the strongest assumptions in our 



approach. It is enough to recall that for a long time the constituent quarks have been considered as a good 
candidate for the correct degrees of freedom in the entire range of energy. The only support for such an 
approach, can be seeen in the success of the Heildelberg group (see Refs. [^] and references therein), in 
the description of soft interactions using this ansatz. 

In the case of mesons, we have only one colorless dipole and we take the transverse wave function in 
the form of a simple Gaussian, namely 

1 

|^'meso„(r)P = -^e (4.12) 

where Sm is a parameter that can be found from the electromagnetic radii, namely, = 0.66 it 0.01/m 
and Rk = 0.58 ± 0.04/m |o|. Using Sm = ^/^Rm we obtain = 1.08/m and Sk = 0.95/m. 

However, electromagnetic form factors gives 
(Q + if = us the space distribution of the electric charge 
inside of the hadron, while in our case the dis- 
tribution of the density of partons (dipoles) is 
Figure 5: Electromagnetic (Fig. ^a) and two gluon (Fig. ^ probed by the two gluon interaction (see Fig. |5|- 
b) form factors, for a meson. Fig. |-c shows the interaction b). This distribution could be different from the 
of a virtual photon with a parton. charged one. We prefer to find the value of Sm 

from the experimental data and, therefore, Sm 
as well as Sp (see below), will be the only fitting parameters in our description of the soft experimental 
data. 

For the baryon, the situation is more complicated: we can have two or even three dipoles. Follow Ref. 



photon 





[23 1, we choose the proton wave function in the simple form for the two dipole model 



lVoton(n,r2)|2 = e ^ (4.13) 



where two dipoles are defined as fi = i?i — R2 and r2 = R3 — ^ \^Ri + i?2 j where Ri is the position of the 
constituent quark i. In Eq. (|1|) Sp = ^J^Rp = 1.05/m for Rp = 0.862 ± 0.012/m@. 

However, for large Nc it is proven that the baryon consists of Nc dipoles Therefore, we consider 
the alternative assumption for '^proton{ri,r2,r3), namely 



— ^ — 



\^proton{ri,r2,r3)\^ = . „ .3 e (4.14) 
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where dipoles have the size Ri — ^{Ri + R2 + Rs)- In this case, from the electromagnetic radius of the 
proton, fohows the value of Sp = Rp = 0.862 it 0.012/m. 

In Eq. ( 4.13| ), we take the size of two dipoles to be equal. We did this for simplicity, since even in the 
constituent quark model (CQM), < |r|| >= 4/3 < \rf\ >. Our hope is that the hadron interaction will be 
determined by the behavior of the dipole amplitude in the saturation domain, where the sensitivity to the 
size of dipoles is expected to be weak. On the other hand, the experimental ratio cjtotij^ — P) / ^tot{p — p) ~ 
2/3 in the CQM stems from quark counting. In our approach, the dipole counting leads to the value of 
this ration 1/2 (for the two dipole model for a proton) or even 1/3 (in the three dipole model) which 
contradicts the experimental data. The only way to obtain a reasonable description, is to hope that the 
perturbative QCD dependence of the dipole cross section on the size of the dipoles a oc will remain in 



the entire accessible range of energies. Fortunately, Ref. |24] demonstrates that this is the case, and this 
result encourages us to search for the description of the soft processes using the dipole hypothesis for the 
hadronic wave functions. 



4.3 Energy variable in the model 

In section 2.2 and 2.3 we have discussed our approach introducing the typical energy variable for deep 
inelastic scattering: xsjorken ^ x = Q'^/s where is the photon virtuality and s is the energy. However, 
with this variable we cannot discuss the soft processes which have = 0. We reconsider the derivation 
of the Bjorken variable to introduce a new energy variable which will have a limit Xsoft — > x at ^ //^ 
where /i is the scale for the soft processes. From Fig. |5|-c we have a relation 

(q + k)'^ = -Q'^ + Xsofts + k"^ = 0; Xsoft = {Q^ + kl) /s (4.15) 

where q is the momentum of virtual photon (q^ — Q"^)- 



In Eq. (4.15) we used the fact that at high energy, /c^ = —k\. Since in the saturation domain 
k\ = Q^(x), we obtain the final expression for our energy variable 

Xsoft — ^ ) Xsoft — *■ XBjorkenj {^■^^) 

S S S 

We believe that one of the main advantages of the saturation approach is the natural choice of the energy 



variable given by Eq. (4.16), which depends on the new scale: the saturation momentum. 



The value of the saturation momentum in our model, we find by resolving the following equation: 

^ {Xsoft, r^sat b = 0) = 1 with r^sat = 4/Ql {Xsoft) (4.17) 

Partons (dipoles) with size rsat are populated densely in the hadron disc, and Qs gives the new scale 
which shows that at < the hard process reaches the saturation domain. 
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Figure 6: Total cross section for proton( antiproton) -proton collision taking into account the interaction of two 
dipoles. 



4.4 Final formulae for proton (antiproton) - proton scattering 



Using the form of the proton wave function given by Eq. (4. IS), we can rewrite Eq. ([4.1D - Eq. ( [4.4[ ) in a 
more accurate form, taking into account the simultaneous interaction of two dipoles (see Eq. (^)). For the 
total cross section, the general formula has the form which follows directly from Fig. & namely 



(Ttot = ^^^^^^1 ^^^2 l^proton. (n,r2) 



(4.18) 



-1 -1 -1 -1 

N{x,ri;b - - f2) + N{x,r2;b - -f2) - N{x,ri;b + -f2) N{x,r2;b + -f2] 
6 b 6 b 



where b is the impact parameter or the distance between ^{Ri + R2 + R3) and the position of the target 
nucleon. One can see that we need to replace 



\'i'hadron{{ri})\'^ N{x,ri;b) 



\'^hadroniri,r2)\' 



(4.19) 



i=l 



-1 -1 -1 -1 

N{x,ri;b - -f2) + N{x,r2;b - -f2) - N{x,ri;b + -r2) N{x,r2;b + -r2] 
6 b 6 b 



in all of the equations of section 3.1. 

Using the explicit form of the proton wave function of Eq. ( [4.13 ), we can easily rewrite the integrals 
in Eq. ( ^ ) - Eq. (|^, in the form 



j (fb(fri(fr2\'^proton{ri,r2) 



(4.20) 



N{x,ri;b -^f2) + N{x,r2;b - l:f2) - N{x,ri;b + ^f2) N{x,r2;b + l:f2) 
6 b 6 b 



drj db\ db\ exp 

Op 



^2 ) ^0 I ^bi 62 

Dp J \Dp 



X {N{x,ri;bi) + N{x,ri;b2) - iV(x, n; 61) A^(x, n; 62)) 
In Eq. ( poD , we used the formula 8.431(3) of Ref. |61]. 
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Since the typical impact parameter increases with energy, (see the next section) and the is smah 
(about 1/3 /m ), we can safely neglect the shift in the definition of the impact parameters in different 
amplitudes, replacing 



|^'Wron({ri})P N{x,ri;b) 

i=l 
1^- 



hadron 



|*hadron(n,r2)P 1 



exp 



-- [Q. [x,ri;b] +9. [x,r2]b 



(4.21) 



(4.22) 



This simplified formula works quite well in all the observables of section 3.1, except for dcTei/dt at large 
values of t. The formula of Eq. ( 4.22| ) has a simple generalization to the case of the three dipole proton 
model . In the three dipole model, for a proton the cross section has the form: 



(4.23) 



(fb(frid^r2 \^proton{ri,r2,n) 



exp 



( X, ri; 6 ) + (x, r2; 6 ) + ( x, ra; 6 



It should be stressed that the second term in Eq. ( 4.18| ) plays a very important role at high energies, 
since it provides the observation that the cross section approaches a black disc limit. Indeed, only with this 
non-linear term in the region where the dipole amplitude tends to unity, (A^ 1) the factor in parenthesis 
(2A^ — A^^) also approaches 1 leading in the black disc regime. 



4.5 Description of the soft cross sections 

In the previous sections we have built our model and fitted all needed parameters using the DIS experimental 
data. In this section, we compare the model with the experimental data without any additional fitting 
parameters. 



4.5.1 Total cross sections 



Using Eq. ( fl.lj ) , Eq. (4.12) and Eq. ( [4.13| ) we calculate the cross sections of pion-proton, kaon -proton. 



proton-proton and antiproton - proton scattering. The results are presented in Fig. |^ and Fig. ^. In 
addition to Eq. ( |4.1| ), we add a contribution of the secondary Regge trajectories using the Donnachie and 
Landshoff parametrization [45|. The agreement with the data is good, and shows that the saturation can 
replace the soft Pomeron, which has been used for fitting the soft experimental data. However, as has 
been mentioned, we still have one fitting parameter, namely, Sm and Sp in the hadron wave function. The 
values of St^^Sk and Sp is given in Table 2. It is interesting to notice that the relation between the fitted 
values of 5", is the same since it stems from the electromagnetic radii. 
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Figure 7: Total cross section for tt — p (Fig. 0-a) and if -p (Fig. |-b) as a function of energy. The solid line (DIS) 
corresponds to the model A which describes the DIS data with -y^ jd.o.J . < 1, while the dotted line presents the 
model B (soft) with a worse x^/d-o.f. ~ 3, but which leads to a better description for DIS, and soft processes data 
together. The number of dipoles in the proton wave function is shown in the legend. 



Model 


Sn (fm) 


Sk (fm) 


5p(fm) 


A(2D) 


0.515 


0.552 


0.561 


A(3D) 


0.515 


0.552 


0.458 


B(2D) 


0.556 


0.597 


0.615 


B(3D) 


0.556 


0.597 


0.519 



One can see that we are able to describe both the value 
and the energy behavior of the total cross section for meson- 
proton and baryon - proton interaction in our two models. 
However, the model B leads to a better description demon- 
strating, in our opinion, that the main idea to replace the 
phenomenological soft Pomeron by the behavior of the 'hard' 



Table 2: Values of parameters m hadronic amplitude in the saturation region, is fruitful . As has been 
wave function that give description of the soft discussed in the previous subsection, the most surprising fact 
data. is that the model describes both the proton-proton and meson- 

proton cross sections. 

This fact can be explained only by the dependence of the dipole cross section, on the size of the dipole 
which the model reproduces. Since such a sensitivity can be only outside the saturation region, or close to 
the boundary between the saturation domain, and the perturbative QCD region. Therefore, we conclude 
that the soft data allow us to obtain new information on this very important transition region between 
the saturation and perturbative QCD region. It should be noticed that we also reproduce the difference 
between the interaction of pions and kaons. Since our QCD dipole cross section does not depend on the 
mass of the quarks, this difference also reflects the difference in the size of the dipoles, inside pion and 
kaon. 

4.5.2 Elastic cross section: energy dependence 

We concentrate our efforts on the proton(antiproton) interactions since meson-proton collisions have been 
studied experimentally only in the limited energy range (see Fig. |8[). The energy dependence of the elastic 
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Figure 8: Total cross section for pp scattering versus energy in our approach. Solid lines correspond to the model 
A, which describes the DIS data with /d.o.f. < 1, while the dotted lines present the model B (soft) with worse 
X^/d.o.f. « 3 but which leads to a better description for DIS, and soft processes data together. 



cross section is presented in Fig. |9|. We use Eq. ([4. 21) for performing our calculations. In Fig. |9| - a we 
compare with the experimental data our model B, while in Fig. ^ - b, is plotted our prediction for model 
A (three quarks). One can see that we obtain a good agreement with the experimental data, without any 
fitting parameter. 

4.5.3 Elastic cross section: t-dependence 



Using Eq. (4^) we calculate the slope for the elastic cross section (see Fig. 10). One can see that in 
the model B, we reproduce the value and energy behavior of the slope at high energy, but overshoot the 
experimental value of this slope at sufficiently low energy. We did not include the secondary Reggeons, 
since we need a piece of information on the slope of the contribution of the secondary Reggeons. We 
also need to take into account the fact that our profile function S{b) corresponds to a power-like form 
factor in the t representation, and the value of the slope for such a function is sensitive on the range 
of t, where the slope was found experimentally. We calculate the slope at t = 0, while the slope at 
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Figure 9: Elastic cross section for proton-proton scattering versus energy in our approach. Solid line shows the 
comparison with the model A with 3 dipoles in the proton wave function which described DIS with ^ jd.o.f . < 1 
while in the dotted line we plot the comparison with the experimental data the model B with /d.o.f. « 3 for DIS 
data for the two dipole model of proton. 



an average value of |t| = |to| > is less than at t = 0. Therefore, it is better to compare with the 
experimental data the t dependence of the differential elastic cross section, given by Eq. (|4.3| ). Such a 
comparison, one can see in Fig. |ll], where the data at small t is described by our model even at low 
energies {y/s = 23.5GeV and G2.5GeV). Notice that at such low energies, the slope, that we predict, 
exceeds the experimental one (see Fig. |lO|). However, in model A, the slope is too small and, actually, this 
is the key difference between the two approaches. The origin for such a difference in the framework of our 
approach is clear: a different energy dependence of the dipole-proton amplitude in the transition region 
between the saturation domain, and the perturbative QCD region. Indeed, Fig. Q shows that model B has 
a steeper behavior in this region in comparison with model A. 



In Fig. IC we compare our prediction with the standard phenomenological parametrization for the slope 



in the soft Pomeron model: B^i = Bq + 2a'ln(s/so) with a' = 0.25 GeF'^ H and Bq = dGeV'^. One 
can see that our prediction in model B is in perfect agreement with this model at high energies, while we 
predict a large slope at low energies. However, model A leads to a smaller slope than is seen in experiment. 
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Figure 10: Slope for the cross section for proton-proton scattering versus energy in our approach. Notations are 
the same as in Fig. ^. 

It should be noticed, the comparison with the experimental data (see Fig. [TT| ) shows that model A fits the 
experimental data. 

In Fig. |ll| we plot the t-dependence for the different energy. We would like to draw your attention to 
two interesting features of our model B: it describes the experimental data quite well at small values of 
t < 0.6 GeV^; and it predicts the minimum at the LHC energy at rather small values of t (\t\ ~ 0.45 GeV"^). 
The appearance of this minimum shows that our model reproduces the effective shrinkage of the diffractive 
peak which results in moving the typical diffractive minima in the region of smaller t, in comparison with 
lower energies. Since we did not include the real part of our scattering amplitude, we see a deep minimum. 
However, the real part of the amplitude will make our minima more shallow. Model A leads to a smaller 
value of Bel, which actually leads to a dip at large values of i, as compared to model B. For our rather 
crude models, such a description looks successful and , we hope, will encourage others to look seriously 
into attempts involving the description of the experimental, data based on the saturation regime of QCD. 



4.5.4 Diffraction production: energy behavior of the cross section 

In Fig. 13 the total cross section (integrated over the masses of the produced particles) of diffractive 



da /dt (mb/GeV^) 

el 




t (GeV^) 



Figure 11: t-dependence of the elastic cross section for different energies a.tt< O.QGeV^. The soUd hne corresponds 
to Model A predictions, while the dotted line describes Model B behavior of the differential cross section. Both models 
reproduce small t behavior of the experimental data in spite of the poor description of the elastic slope in model A. 



production is plotted. We use Eq. ( |4.6| ) to calculate this cross section. However, we multiply this formula 
by 2 since the two protons: projectile and target, can dissociate diffractively. One can see from this picture 
that we reproduce values and energy dependence in agreement with the experimental data. The qualitative 
behavior of the experimental data on Cdiff versus energy is characterized by the slow dependence on energy 
in the region of high energies. This fact is perfectly reproduced by our model. 



However, one can see from Fig. 13 that our approach cannot reproduce asD in the region of low energies 



(small values of produced mass). There are at least two reasons for such a failure: first, we do not take into 
account the exchange of the secondary Reggeons whose contribution is essential, as we see on the example 
of the energy behavior of the total cross section, and, second, the contribution of the one extra gluon 
emission has been calculated in leading log(l/x) order, which cannot describe the low mass diffraction. 

|]) can be described assuming, in addition to 



The experimental data on single diffraction (see Refs. |53 



the triple Pomeron vertex, that we modeled by extra gluon emission, a significant contribution from the 
Reggeon-Pomeron-Reggeon vertex in Regge phenomenology, which we cannot include in our approach. 
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Figure 12: Total cross section for meson-proton and proton(antiproton)-proton scattering scattering at high energies 
(extrapolation). Notations are the same as in Fig. ^ 



4.6 Predictions for the LHC range of energies 

As we have seen, the model is able to describe the available experimental data quite satisfactorily (and 
Model B even quite well) both for soft and hard interactions. Therefore, we can rely on the model B for 
the predictions in the LHC energy range. From Fig. ^Fig. 13 we see that at the LHC energy we have 
atot = 101.3 m6, aei = 28.84 m6 , Bei = 18.4 Ge^-^ and (Xdi// = 10.5 m6 . 

Table 3 shows that our prediction for LHC energy is close to ones that are given by the models that 
fitted the experimental data. The RMK and GLM2 models are based on soft Pomeron phenomenology 
while the GLMl model is close to our approach: in this model the existence of soft Pomeron is not assumed 
and the parametrization was chosen for from the same ideas as our approach here. However, in the GLMl 
model all parameters are found from fitting soft data. In spite the fact that the values of cross section are 
close the different models are different in more detailed characteristics. For example, one can see from the 
Table 1 they predict different values for Af,i{b = 0) . From unitarity Aei{b = 0) < 1. Our model predicts 
that at the LHC energy our proton-proton collision is close to the black disc regime. The same we can see 
in GLM2 model, but in the RMK model we are not so close to this regime at the LHC, while in the GLMl 
models the proton-proton interaction is far away from the black disc limit. 
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Figure 13: Cross section for single diffraction production for proton(antiproton)-proton scattering versus energy. 
The upper curves are the total cross section for diffractive production while the low curves show the energy behavior 
of the low mass contribution to the diffractive cross section from Eq. (4.5). Notations are the same as in Fig. |9[ 



Model 




0"ei 




Bel 


Aei{b = 0) 


Aei{b = 0) 




mb 


mb 


mb 


GeV-^ 


(Tevatron) 


(LHC) 


Our A(3D) 


83.0 


23.54 


10 


16.67 


0.94 


0.98 


Our B(2D) 


101.3 


28.84 


10.5 


18.4 


0.94 


0.98 


GLM1[||] 


110.5 


25.3 


11.6 


20.5 


0.6 


0.7 


GLM2[||] 


91.7 


20.9 


11.8 


17.3 


0.94 


0.95 


RMK|| 


88.0 (86.3) 


20.1 (18.1) 


13.3 (16.1) 


19 


0.89 


0.92 



Table 3: Cross sections and elastic slope at the LHC energy in different models. For the RMK model, we put two 
parameterizations and the value for B^i directly from the curve in Fig. 18 of Ref B3| 



As you can see from Eq. (§)-c we have some information on the total cross section from the cosmic ray 
experiment. However, to extract the value of the total cross section from such an experiment, we need to 
know the meson-proton cross section as well. In Fig. we plotted our predictions which we hope will be 
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useful for discussing the cosmic ray experiment*. 

5. Survival probability for diffractive Higgs production 




Figure 14: Diffractive Higgs production. 



The diffractive Higgs production, is the reaction which has the best experimental signature for the 
discovery of the Higgs boson at the LHC. At fist sight this process occurs at short distances of the order 
of I/Mh, where Mh is the mass of Higgs boson (M^ > lOOGeF) and can be calculated in perturbative 
QCD (see Fig. 14-a and Ref.|46| for calculations). However, as was noticed a long ago E^, El] that it is not 



enough to calculate the diagram of Fig. 14-a which describes the Higgs boson production from one parton 



shower. We have to multiply this cross section by a probability that two parton showers (or more) will not 
interact with the target since they could produced hadrons that will fill up the large rapidity gaps between 
Higgs and protons in the final state (see Fig. |l^b). This probability we call the survival probability. In 
the case of our model with the eikonal type formula for the scattering amplitude, the expression for the 
survival probability is very simple (see Fig. [ij-c), namely, [47, 48 



(l^'l) 



Jd% (/tiVi(i2r2|^p™ton(ri,r2)|2e~^^(^''^i'''))' anib-Fig. 
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(5.1) 



The appearance of the factor exp (— ^r2(s, ri, 6)) in Eq. (5.1) is clear from Fig. |T^c) but the power of 4 
requires discussion. Actually this power reflects the fact that we have to find the probability that neither 
dipole 1 nor dipole 2 could scatter inelastically. The probability that one dipole does not scatter inelastically 



*We thank S. Nussinov for drawing our attention to tlie necessary knowiedge of tfie meson-proton total cross sections at 
high energy. 
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is equal to 



e 2 



Therefore for the probabihty that two dipoles cannot scatter inelastically we obtain ({\e 2^('*'^i''')|)^ . We 
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Figure 15: The value of the survival probability for diffractive Higgs production versus energy. Solid lines show 
the survival probability in Model A that describes the DIS data with a very good /d.o.f < 1, while dashed lines 
correspond to < S*^ > in Model B which gives a good description of all soft data. 



can see from Eq. (|5.1|) that for the calculation of the survival probability we need to know the b - dependence 



of the hard cross section. This dependence has been discussed in detail in Ref.|49], where it was extracted 
from the process of diffraction production of J/^' with the proton (elastic) and the state with mass larger 
than the proton (inelastic). These two processes have two different slopes in the t behavior: B^i = AGeV^'^ 
and Bin = 1-86 GeV"'^^^ . The dependence of the vertex for photon to ,1/^ transition was extracted from 
these reactions in Ref. leading to the values of slopes for the vertices of the transition of proton to 
proton (elastic), and proton to non-proton final state (inelastic): Bg.i = 3.6GeV~^ and Bin = lA6GeV~^. 
In our model, the projectile proton first goes to the state of two free dipoles (see Fig. |l^c) , these two 
dipoles can produce two and more parton showers; and finally two free dipoles in hard processes create a 
proton. Therefore, at least the upper vertex is the same as in the inelastic diffractive production of J/\I' 
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in DIS. Considering the lower vertex being elastic, we can write cr^ib; Fig. 14 — a) in the form 



aH{b;Fig. 14 



ttRh 



(5.2) 



with = Bin + B(,i = 5.06 GeV^^. However, it is our model disadvantage that we describe differently 
projectile and target protons. In any case the eikonal type model that we use, does not mean that only the 
elastic rescattering contributes to the shadowing corrections. Treating both protons on the same ground 
we take R% = 2Si„ = 2.92GeV^'^. 

Since our model naturally includes both soft and hard interactions we can hope that the hard process 
of diffractive Higgs production can be described in the same way as we did in the model using the same 
amplitude O. However, we need to take into account that in this case 



aH{b;Fig. 14 



aoS\b) 



(5.3) 



where S{b) is given by Eq. ( 3.15| ) 



Our results on the survival probability is shown in Fig. |15|, and in Table 2. Table 2 demonstrates that 
our value for the survival probability turns out to be much smaller than the values in the phenomenological 



models, that fitted the experimental data using the soft Pomeron approach (see Table 2, the review of |54] 
and Ref. |^5|). We think that the difference stems not from the details of the models, but from the key 
ingredient of our model: the fact that we took into account the interaction at short distances. Therefore, 
we confirm that the short distances give a substantial contribution to the value of survival probability, 
as has been noticed in Refs. 58]. The Table 2 also shows that the value of the survival probability 
crucially depends on the model for the 6-dependence of the hard cross section. Our model with Eq. (l5 



for the hard cross section is similar to the one that was used in the GLM model |25]. It is interesting that 
our model confirms the general tendency to obtain a smaller value for the survival probability advocated 
in Ref. Mi. 



6. Lessons from the model 



Long distance physics is very complicated, non-perturbative phenomenon, and we certainly do not pretend 
that we are able to describe it in its full richness. However, we demonstrate in this paper that the gap 
between this physics and the short distance physics which is under full control of perturbative QCD, is not 
so huge that it would be a hopeless task to build a bridge. Comparing this with the experimental data, we 
showed that the soft data depends on the transition region between saturation domain and perturbative 
QCD region and, therefore, can give valuable information on this transition, checking our theoretical 
approaches to it. The widely used, phenomenological soft Pomeron do not appear in our approach, and 
we hope that the reader will ask the question: do we need a soft Pomeron, having in mind our negative 
answer. 

Our description of the experimental data, is not worse (in the case of Model B) than the one in the 
models which fitted the data on the basis of the soft Pomeron phenomenology (see Refs. 
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Model 


CJHib) 


(l^'l) 


Our Model A(3D) 

Our Model B(2D) 

GLM1|25| 
GLM2|56| 
RMK[55| 


Eq. Q Rjj = 5.06 GeV-^ 
Eq. Q Rjj = 2.92 GeV'^ 

Eq. dU) 
Eq. dU) Rjj = bMGeV-'^ 
Eq. dU) = 2.92 Gey^^ 

Eq. dU) 
Two channel model for an 
Two channel model for an 
Eq. Q = llGeV"^ 
Eq. Q = 8 


0.24% 

0.02% 

0.89% 

0.24% 

0.096% 

0.57% 
2% (0.7%) 

0.21% 
3.2% (2.3%) 
1.7% (1.2%) 



Table 4: The survival probability for diffractive Higgs production at the LHC energy, in different models. In the 
GLM model, for <JH{b) is used in the two channel model, with the same Bin and i?ei for the diffractive production 
of J/^- in DIS. 



Recalling that we fitted all parameters of our model from DIS processes, we interpret this success in the 
way that the high energy soft scattering processes are determined by QCD, at short distances of the order 
of 1/Qs, where Qg is the saturation momentum. Model A describes the data worse than Model B, but we 
consider this description quite satisfactory remembering the crude character of this approach. It should be 
stressed, once more that the difference between our model A and model B, is in our attitude to the simple 
formula of Eq. ( |3.17| ): in model A we trust Eq. ( |3.17 ) in the entire kinematic region of accessible distances, 
while in model B we view this formula as a kind of qualitative description, that includes the main features 
of the saturation regime. Therefore, we were searching for the parameters of model B, in the way that 
describes all the data both on DIS, and on soft interactions, in the best possible way. 

In our model, we used several assumptions which have a different theoretical status. The first assump- 
tion is the exponential form of the dipole scattering amplitude (see Eq. ( p. 11 ) and Eq. ( |3.17| ). We have 
discussed the theoretical arguments for such an assumption, namely, this form has been proven for the 
transition region between perturbative QCD, and the saturation domain (see also Ref. 10]). It should be 
stressed that the soft data, are sensitive to the transition region as we have discussed. 

The second assumption, is the expression of Eq. ( 3.1(j ) for (x, r, h). This assumption is a compromise 
between what we should do, and what we can do. We can check it by describing the DIS experimental 
data, at large values of the photon virtualities, where the difference between our formula for fi, and the 
DGLAP expression for it should be large. 

The third assumption is the impact parameter dependence oi^l{x, r, b) (see Eq. ( [3.141) and Eq. ( 3.15 )). 
This is a pure phenomenological ansatz, since the current stage of our theory does not allow us to find 
b dependence |^. Actually, only soft data allows us to check this dependence. Indeed, our $7 describes 
the t-dependence of the elastic cross section while , for example, the Gaussian b dependence results in the 
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appearance of the structure of maxima and minima, at small values of t, which contradicts the experimental 
behavior of duei/dt. 

In Fig. [l6[-a we plot the dipole amplitude averaged with the parton wave function 



(|iV|)l 



(fri(fr2\'l!proton{ri,r2)\^N{x,ri,b = 0) 



(6.1) 



One can see, that this amplitude increases and it approaches 1. However, it happens at ultra high 
energies and at the Tevatron energy, for example, this amplitude is only 0.8 at 6 = 0. Such an average 
scattering amplitude, leads to the elastic amplitude of proton-proton scattering 2N — N'^ = 0.96. The 
averaged slowly increases with energy, and reaches 1 at energy. At the LHC we expect (|A^|) = 0.915 and 
the average elastic amplitude is close to unity (|2A^ — A^^|) = 0.99. Such a behavior shows, that the values 
as well as the dependence on energy crucially depends on the behavior of our dipole amplitude, in the 
vicinity of the saturation scale. To illustrate the strength of the saturation effect, we plot the average $7 
at 6 = as a function of energy (see Fig. 16 -b), which is a considerably overshoots the average scattering 
amplitude. It should be stressed, that in spite of the fact that the elastic amplitude is very close to unity 



at the LHC energy, one can see from Fig. 16-a that (|A^|) is only 0.9. and the asymptotic behavior with 
(|A^|) close to 1, starts from s = 10^^ GeV"^. 




\ogJ s/1 GeVO 



Fig. |16|-a 



\ogJ s1/GeV0 
Fig. H-b 



Figure 16: Average dipole amplitude {{N) see Eq. ( |6.l[ )), average proton- proton amplitude (both in Fig 
average opacity r2(see Fig. |l^b) at 6 = as a function of energy in our model. 



13-a) and 



The principle difference between our models and the models of soft interactions, is the fact that 



we predict all observables fitting all parameters from the DIS data. The GLMl model|25] is close to our 
approach ideologically, because it does not assume the existence of the soft Pomeron. However, at first sight 
this model has two major shortcomings: minima at small t which have not been seen experimentally; and 
the slow fall down of the cross section of diffractive production, as a function of the mass of produced system 
of hadrons. Our model shows, that the t dependence can be easily heeled by assuming the exponential 
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form for the profile function, instead of the Gaussian one that has been used in the GLMl model. As far 
as large mass diffraction that has been neglected in the GLMl model|25|, we show that this diffraction is 
essential. In this respect, we are close to the RMK model and to the GLM2 model [56| which are based 
on the soft Pomeron exchange. Our model can be considered as an argument that the multi Pomeron 
exchanges, and the Pomeron interactions could be essential for high mass diffraction. 

In general we demonstrated in this paper that the distances, essential in so called soft interactions, 
is not so long, but rather about l/Qg Rh, where Rh is the hadron radius. This conclusion stems both 
from the success of our model in the description of the data, using the parameters fitted in DIS, and from 
the use of the energy variable Xsoft^ which is determined by the saturation scale. 

We hope that our model will generate deeper theoretical ideas on the matching between soft and hard 
interactions, based on high parton density QCD. 
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